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We review an evolution equation for quantum entanglement for 2x2 dimensional quantum systems, the 
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I. INTRODUCTION 

In quantum mechanics, the dynamics is introduced by 
virtue of the Schrodinger equation: 



ft*W»=i*W)>; 



(i) 



where |v(0) represents the state of the system at time f, and 
H is the Hamiltonian governing the evolution. For a time- 
independent Hamiltonian, solving the dynamics of a quantum 
system "merely" requires to diagonalize its Hamiltonian in or- 
der to arrive at the time evolution of a set of basis states - the 
eigenbasis. Since the Schrodinger equation is linear with re- 
spect to the state, this immediately also provides the time evo- 
lution of any initial state, by expressing the latter in the eigen- 
basis of the Hamiltonian. Thus the knowledge of the Hamil- 
tonian's eigenstates suffices to construct the evolution of any 
other state as well as for every quantity which is a function 
of the state, e.g. fidelities, overlaps, and expectation values of 
observables. However, for large quantum systems, diagonal- 
ization of the Hamiltonian quickly turns into a tedious ven- 
ture. It is thus mandatory to directly construct an evolution 
equation for the quantities of interest, which comprises only 
the minimal set of information needed. 

Entanglement, arguably the most pronounced signature of 
quantum mechanics, is an even more extreme example of the 
intractability mentioned above. Entanglement is a non-linear 
function of the state from its very definition. A state p of a 
quantum system that is identified with the Hilbert space J#a <8> 
,3%b is said to be entangled, if and only if it cannot be written 
as a separable state, i.e. in the form (T): 



(2) 



where pi > and p, ■= 1. Therefore, in order to quantify the 
amount of entanglement inscribed in a state, one must search 
for a decomposition of p which is the "closest" to a state with 
the above form - where the definition of "close" depends on 
the entanglement measure being used. One possibility is the 
convex-roof construction |2 |, by which the entanglement mea- 
sure E defined over the pure states is extended to mixed states 

P=LiPi\Wi){Wi\^T- 



E(p) = inf 2>£(|W»; 

{Pi,\¥i)}j 



(3) 



where the infimum is performed over all possible decompo- 
sitions of p (see Ref. [3| for other possible definitions). Two 
important characteristics of this expression must be noticed: 
i) E is in general a non-linear function of the state since it is 
convex: E(Xpi + (1 - A)p 2 ) < lE(pi) + (1 - k)E(p 2 ), and 
if) it is hard to be evaluated due to the infinitely many decom- 
positions. In fact, it was shown in |4] that in general to decide 
if a given state is entangled or not is a NP-hard problem. 

The difficulty to determine the amount of entanglement is 
magnified when the quantum system undergoes some dynam- 
ics. Up to recently ]5]|6) this required to follow the evolution 
of the state and subsequently, from the knowledge of all its 
parameters, calculate its entanglement (see Fig. ITJ. Experi- 
mentally this implies a full state tomography for each point 
in time, and for each initial condition. This is specially de- 
manding in realistic situations where quantum systems inter- 
act with uncontrolled or unobserved degrees of freedom, usu- 
ally lumped together under the term environment. The dy- 
namics due to a coupling to an environment usually degrades 
the quantum features of the considered quantum system, and 
it is often exploited to explain the classical behaviour of the 
world despite it being quantum on the atomic level [7|. Fore- 
most, in situations where true quantum effects such as entan- 
glement provide a resource, e.g. for teleportation and various 
other quantum information and communication tasks O, it is 
of primary interest to know how such resource decays while 
being processed or stored. 

On the following pages we develop a deterministic equation 
for the evolution of quantum entanglement. This relation con- 
tains only a single quantity which, independently of the initial 
state, describes how a state's entanglement evolves. This sin- 
gle quantity constitutes a benchmark for given open system 
and decohere nee dynamics. 

We start in section[ll]by reviewing some basic concepts and 
formulae concerning the entanglement measure of our choice. 
After that, in section [Hi] we present an equation of motion 
for the entanglement of an initially pure state of two two-level 
systems when either one of them undergoes some evolution 
(as described by a completely positive map); extensions of this 
setup are also introduced. Finally, in section [TV] we comment 
on the generalization of our framework to higher dimensional 
systems. 
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E[p(0)] 



P(0) 



entanglement evolution 



state evolution 



-E[p(t)] 



p(0 



where A,- are the eigenvalues of the matrix p (<7 V <8> °y)p* (°y ® 
<jy) with A] denoting the largest among them, and the complex 
conjugation is taken with respect to the product basis of eigen- 
vectors of a z . Although this analytical expression for mixed 
states enables and encourages the treatment of 2 x 2 systems 
under exemplary cases of open system dynamics, we rather 
focus on the algebraic property introduced in ([5]) to derive our 
results and build intuition. 



B. Entanglement under filtering operations 



FIG. 1: From state evolution to entanglement dynamics. A direct 
equation of motion for entanglement, which uses only the minimal 
necessary information about the dynamics and initial state is desir- 
able. 



II. BASICS 



In preparation of what follows, we briefly review how en- 
tanglement changes under filering operations lfl3l fT4l . Such 
operations can be understood as the result of a measurement 
in a higher dimensional space. We consider a quantum system 
composed of two subsystems with state \yr) = £i jWij\i)\J) 
and apply a filtering operation M = Y,ki Mia \k)(l\to one of the 
subsystems. After renormalization the state becomes 



Concurrence 



We choose concurrence as a quantifier of entanglement Q 
because of its convenient algebraic properties. For a pure state 
\x) of two finite-dimensional quantum systems, say with re- 
spective Hilbert spaces M'a and Mb, it is 



C{\%)) 



1-Trp 



A/B 



(4) 



where p A / B = Tr B / A \x)(x\ is me reduced state of either sub- 
system after tracing out the other. The compound quan- 
tum system of the smallest possible dimension which can ex- 
hibit entanglement is composed of two subsystems with two 
levels each. The corresponding pure state, when expressed 
in an arbitrary product basis of the respective subsystems, 
\x) = L, 1 j=oXij\ e i) \fj)i yields the expression for concurrence 
by means of a determinant of its coefficient matrix %\ 



C(|z» = 2|det*|. 



(5) 



When the quantum system undergoes open system dynamics 
and couples to its environment, the quantum state becomes 
mixed in general. The resulting mixed state p is described by 
a probabilistic mixture of an ensemble of pure states, i.e., p = 
HiPi\Wi)(Wi\ w i tn corresponding probabilities /?, that sum to 
one. As previously stated, this measure can be extended over 
the mixed states by virtue of the convex roof construction |2 
ED: 



C (p) = f inf I>C(|W>)- 

{Pi,\¥i)}i 



(6) 



In particular for states of high dimensional systems this opti- 
mization is difficult to perform. In fact, analytical expressions 
are only known for states of high symmetry ifTTl [T2l and for 
systems of dimension 2x2. In the latter case, the minimum 
is obtained by [ 9 1 



C(P) 



{O, 



(7) 



IV, 



M(g>I\y/) 



: Ev(#>u>, 



(8) 



where the new matrix elements with respect to the product 
basis are yr^. Y^i^kiWij' apart from the normalization factor, 
i.e. the product of the matrices: yf' °< M\j/. Since for a ma- 
trix product the determinant fulfills det(My/) = detMdety/, 
the concurrence |5} of the pure state after a filtering operation 
changes to 



„ , 1N IdetMl „ . , 

C ^ = llJo Tl M|2 C y • 

||M<x>l|y/)||- 



(9) 



If M is invertible, e.g. it is not a projective measurement, the 
entanglement is not completely destroyed but merely rescaled 
by the determinant of M. Thus, for suitable combinations of 
filtering operations M and states | yf) entanglement can proba- 
bilistically increase. 

Now, the treatment for pure states allows us to extend the 
above relation to mixed states. A filtered mixed state is 



P = 



(M<g>I)p(M<x>II) + _ (>f<g>ll)p 



Tr[(M<g>I)p(M<g>]I) + ] Tr[(>f <g>I)p] 



(10) 



where we abbreviated the action of M onto p as the action of 
a map Therefore, the pure states of the decompositions of 
p and its "filtered" version p' are only related by the filtering 
operation. This means that if the infimum in |6| for p is at- 
tained for a mixture of pure states { | y/j) } with respective prob- 
abilities pj, then the optimal decomposition for p' is realized 
by {ly/)} which are the filtered versions of the with 
the same probabilities. Therefore we obtain for the change of 
concurrence of a mixed state under a filtering operation 1 1 5 1 : 



C(P') = 



IdetMl 



Tr[(^®I)p] 



C(p). 



(11) 
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HI. EVOLUTION EQUATIONS FOR CONCURRENCE 

In this work we model the effect of open system dynam- 
ics and in particular of decoherence on the quantum state by 
means of dynamical maps. This provides the required free- 
dom to treat general open system dynamics without having to 
assume a specific interaction. Such maps act on density oper- 
ators and mediate the evolution that the state undergoes during 
a certain time interval ll8l [T6l : 



p(0=A[p(0)]. 



(12) 



The map A is linear, preserves the trace, and maps quantum 
states to quantum states. In order to describe the dynamics 
of systems which are part of a larger system correctly, A also 
needs to be completely positive. Any such map can be written 
as p(0) h- > Li ^i P(O)^/, where K( are called Kraus operators 
and, due to the trace-preserving constraint, fulfill -kJK{ = I. 
This dynamical map is usually termed Kraus or operator-sum 
representation (SJ. 

We begin the treatment of the concurrence evolution under 
such dynamics with the simplest possible case, i.e. a 2 x 2 
system that starts in an initially pure state \%) and of which 
only one subsystem undergoes open system dynamics. A sce- 
nario that approximates this setup is a laboratory that prepares 
these two quantum systems in an entangled state, and subse- 
quently sends one of the systems to a remote laboratory, for 
instance through a glass fiber in case of photons. During this 
process the transmitted photon suffers the detrimental influ- 
ence of its environment, designated by A, whereas the other 
quantum system remains well protected in the lab. The sys- 
tem's final state thereafter is given by 



p = (I®A)\x)(xl 



(13) 



We aim for a relation between the final and initial amount of 
entanglement, C(p) and C(\x)), respectively, in terms of A. 

In order to derive such an evolution equation, we reformu- 
late the initial setup ( p~3] > equivalently into its dual setup, i.e. 
we want to exchange the role of states and channels. In the 
first step, we express the initial state \%) = jXij\ ') \j) as me 
result of a filtering operation M onto the first subsystem of a 
maximally entangled state: 



\X) =M®I|0- 



(14) 



where we choose \<j) + ) — (|0)|0) + |l)|l))/\/2 as an ex- 
emplary maximally entangled state. We thus define 
M = Vi^k i Xkl\k) {l\, which already includes the normaliza- 
tion factor. Therefore, the same final state of the original setup 
( fT3| l can also be prepared via A)|0 + )(0 + |. Since the 

filtering J£ and the open system dynamics A act on differ- 
ent subsystems, their order may be interchanged. Applying A 
first, we arrive at the dual setup: 



p = {Jl <8>I)Pa, 



(15) 



where Pa = (I<8> A)|0 + )(0 + | is a mixed state in general and 
contains all the parameters of the dynamics A. We have thus 



interchanged the roles of initial state and dynamics by moving 
all the parameters of the initial state \x) into a map and 
the parameters of the dynamics A into the state Pa- This is the 
essence of the Choi-Jamiolkowski isomorphism |[T7l[T8l from 
which the state Pa is shown to be isomorphic to the map A. 

We exploit this equivalence and apply the filering result 
( fTTj ) to the dual setup ( fT3j ). The resulting determinant of the 
filtering M equals the concurrence of the initial state \x) and 
yields the first result: 

C[(I®A)|Z)(Z|] =C(\x))C[(I®A)\<l> + ){<l> + \]. (16) 

The concurrence of all pure states evolves exactly as that of 
a maximally entangled state, and is merely rescaled by the 
initial amount of entanglement. This relation clearly sepa- 
rates the initial condition from the evolution, which may be 
any dynamical process. The maximally entangled state thus 
serves as a benchmark according to which the entanglement 
of all other pure states evolve. This means that for a given 
open system dynamics A, there is qualitatively only one type 
of entanglement evolution. Also note that the particular kind 
of the maximally entangled state does not matter, since any 
other maximally entangled state, e.g. a singlet state, may be 
prepared from by a single unitary operation on the "left" 
subsystem - which changes Pa but leaves its entanglement 
invariant. The relation holds also for more general A than 
we initially assumed. Maps which do not perserve the trace, 
ZiKjKt < I, also fulfill ([T6) but with concurrence evaluated 
for the resulting unnormalized states. Alternatively, when us- 
ing the normalized versions, the right side needs to be rescaled 
by the respective probabilities of \x) and |0 + ) to pass A: 



C 



Tr[(I® A) 



xC 



C(\X))^X 
Px 

(i®A)14>+)(4>+ 

Tr[(I<g>A)|0+)<0- 



, (17) 



where p x = Tr[(I® A)|^)(^|] and similarly for p§ with |0 + ) 
instead of \x). 

Since pure states are an idealization, let us now consider the 
case of two two-level systems in an initially mixed state po. 
After one of the systems is subjected to non-trivial dynamics 
A, the state reads 



p = (I<g>A)po 



(18) 



In order to obtain a relation similar to ( fTB"] ), but for an ini- 
tially mixed state, we relate the latter with a decomposition 
into pure states and then employ ( |16) , For this purpose we 
use an optimal decomposition of po, namely one for which 
C(po) = Hi PiC(\xd)- Since the map A is linear, the mixture 
of the \Xi) individually undergoing the map constitutes the fi- 
nal state p. However, concurrence is a convex function as 
every valid entanglement measure, meaning that probabilistic 
mixing of states cannot increase the entanglement on average, 
formally 



C 



J>(I®A)|js)<fc| <^iC[(I®A)|^>0t/|]. (19) 
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Therefore, we obtain an upper bound when applying (T6) to 
each of the terms on the inequality's right hand side. After 
collecting the sum over C(\%,)) into the initially mixed state's 
concurrence, we arrive at the evolution equation for mixed 
state entanglement 

C[(I®A)po] <C(po)C[(I<g>A)|0 + >(0+|]. (20) 

With this relation we can straightforwardly approach more 
elaborated setups where both subsystems undergo arbitrary 
open system dynamics: 

p = (Ai®A 2 )p . (21) 

We group the effective result of just one of the maps onto po 
into an intermediate state, say p = (1(g) A 2 )pi with p\ = (Ai (g) 
I)po, and apply the relation for mixed states ( |20] > twice to yield 

C[(Ai ® A 2 )p ] < C(p ) C[(I(gA 1 )|0+)((/.+ |] x 

xC[(A 2 ®I)|0+)(0+|]. (22) 

This relation describes the evolution of entanglement as both 
subsystems undergo arbitrary dynamics for general initial 
states. Once more, the behaviour of the maximally entangled 
state undergoing the local dynamics separately determines the 
upper bound. Also in this more sophisticated setup the initial 
condition C(po) and the dynamics completely separate such 
that statements about the dynamics, here the evolution of the 
maximal possible entanglement content in the state, do not 
require the knowledge of the initial state. 

Despite the relation being an inequality in general, the 
equality holds in many relevant cases, for example when one 
subsystem of an initially pure state undergoes a series of con- 
catenated dissipative dynamics. 

Relations ([16) and |22) fostered experimental work using 
linear optics with entangled photons from spontaneous para- 
metric down-conversion [22, 23 j. In contrast to upper bounds 
( pO) and ( f22) , when considering entanglement as a resource, 
it is more relevant to know how much of it is at least avail- 
able, i.e. guaranteed to be left after a corrupting environment 
interaction. Such lower bounds are presented by several pro- 
posals lf!9ll20ll2TH. 

Although relations ( pO) and ( p2| ) provide strictly speak- 
ing only upper bounds, we demonstrate here that neverthe- 
less they provide a good approximation for experimentally 
relevant domains. Experiments that aim to create, manipu- 
late, and possibly utilize entangled states usually operate with 
states of a high purity |24|. Also, the time scales of these ex- 
periments fit the time domains where entanglement is mostly 
preserved. In order to profile decoherence effects during rel- 
evant, i.e. short time scales, we only need to focus on ini- 
tially pure states and slight, but non-vanishing, decoherence 
effects. Since all ingredients in the derivation of ( p2) are con- 
tinuous with respect to small changes in the initial quantum 
state, equality in ( p2) holds in the limit of pure initial states 
and vanishing effect of at least one of the maps Ai/ 2 . For 
small amounts of decoherence acting on both subsystems in- 
dividually, Fig. [2] shows (for several examples of maps) that 
in particular for highly entangled states the upper bound gives 



a very good approximation. Since the obtained relations for 
the evolution of concurrence are independent of the particular 
initial state, also a numerical sampling must avoid any bias. 
We thus sample uniformly over pure states l25l . Under such 
uniform distribution the point density in Fig. [2]indicates that 
most states have a decay similar to the bound. This is quanti- 
fied by the histograms in Fig. [3] which show that for a weak 
environment influence the great majority of all uniformly sam- 
pled states evolves very close to the bound. The upper bound 
thus leads to a good approximation for most pure states in the 
considered situations. 



IV. ENTANGLEMENT EVOLUTION OF HIGH 
DIMENSIONAL SYSTEMS 

An extension of the relations ([16) and ( |22) to quantum sys- 
tems of higher dimension is possible [6 1. However, the arising 
complexity in higher dimensions demands its toll. Although 
concurrence |4) is available for systems of say d x d dimen- 
sions, it can no longer be calculated via the determinant of 
the state's coefficient matrix (|5}. An entanglement quantifier 
that can be evaluated for mixed states via such a determinant 
is G-concurrence [26], which reduces to concurrence for 2x2 
dimensions. In particular, G-concurrence quantifies the entan- 
glement which requires all the systems dimensions, and thus 
vanishes for states that live on a strict subspace. For example, 
a Bell-like state (|1) 1 1) + |2) \2))/V2 in a 3 x 3 system gives a 
vanishing G-concurrence despite being entangled. However, 
among all pure states in d x d dimensions, those with van- 
ishing G-concurrence are of volume zero. For a pure state 
\X) = T.tj=i XijV) \J) G-concurrence is given by 

C d (\x))=d\det X \ 2/d (23) 

and it is extended to mixed states by a similar convex-roof 
optimization as concurrence in (|6). Unfortunately, despite the 
simple algebraic form for pure states, an analytical expression 
for mixed states is so far not available. A similar approach 
as the one for concurrence, however, yields upper and lower 
bounds 1271 . 

Since the argumentation for the evolution equations for en- 
tanglement of a 2 x 2 system does not depend on the dimen- 
sion but rather on the determinant structure of concurrence, 
we can straightforwardly apply the same reasoning for G- 
concurrence and obtain an evolution equation of the very same 
structure: 

C d [(l®A)\ X )(x\] =C d (\x))C d [(l®A)\(l> d )(<l> d \], (24) 

where the maximally entangled state in d X d dimensions is 
\<j> d ) = Yf!=i \i)\i)/"vd- When both subsystems undergo open 
system dynamics and/or the initial state is mixed an upper 
bound results 

Q[(Ai ®A 2 )p ] < Q(p ) AOIfcXtol] x 

xC d [(A 2 ®I)\h)(to\]. (25) 
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i.or 






C(p) [upper bound] C(p) [upper bound] 



FIG. 2: Bound tightness for different initial entanglement. 10000 initially pure states of a 2 x 2 system are uniformly sampled, and each 
subsystem evolves under identical incoherent dynamics. The left panels show the data for a phase damping environment in which the subsys- 
tems do not exchange excitations with the environment, but lose the coherences in the respective subsystems with probability p due to elastic 
scattering. The corresponding Kraus operators are Kq = \/\ — pi, K\ = •/p\g){g\, and K2 = y r p\e){e\, where |g) and \e) denote ground and 
excited state, respectively. The right panels show data for a dissipative map where p is the probability for an excitation to be released. The 
corresponding Kraus operators are Kq = \g) {g\ + ^/l — p\e) {e\ and Ky = ^fp\g) (e\ . Points indicate the concurrence of the final state C(p) in 
units of the upper bound (which is tight for these examples of dynamics) versus the initial amount of concurrence C(\%)). Darker regions 
contain more datapoints. 



Although concurrence and G-concurrence can assume values 
largely independent of one another, concurrence cannot be- 
come arbitrarily large without entanglement that involves all 
levels of both subsystems and hence G-concurrence becoming 
non-vanishing. With a vanishing G-concurrence, concurrence 
at most attains ^/l — \/(d— l) 2 of its maximal value. There- 
fore, if concurrence exceeds this fraction of its maximal value, 
G-concurrence is necessarily present. 



V. CONCLUSION 

Despite of its hard characterization, here we explored the 
possibility to design deterministic equations of motion for en- 



tanglement inscribed into two finite dimensional systems. For 
the case where the initial state is pure, and only one of the 
subsystems undergoes an incoherent process, the equation of 
motion for entanglement assumes the form of a simple fac- 
torization law - the first term contains only information about 
how the entanglement is affected by the dynamics, and the 
second term scales the first by the initial amount of entan- 
glement. For more realistic scenarios, where both parts are 
influenced by an environment and/or the initial state is not 
pure, the intricate nature of entanglement only allows for the 
derivation of generic upper bounds. Nevertheless, for most of 
the situations of experimental interest, the bound represents a 
fiducial description of the entanglement dynamics. An appar- 
ent reduction of complexity in determining the entanglement 
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C(p) [upper bound] C(p) [upper bound] 

FIG. 3: Histograms of the final states' concurrence after a phase damping (left) and dissipation (right) map on both subsystems, respectively. 
The concurrence after the evolution is plotted in units of the upper bound. All histograms are normalized such that they approximate the 
probability density and correspond to summing events in Fig. [2] along columns. The initial states are uniformly sampled. The sample size is 
10000 states. 



of a state is obtained when the dynamical process is taken into 
account. 

It is important to note that, the results here presented are 
rooted in the Choi-Jamioikowski isomorphism ifTTl [T8l . and 
on the mathematical characteristics of the entanglement mea- 
sures used. For a more general framework, where the dynam- 
ics are not restricted to local ones, the entanglement measures 
do not present favorable mathematical structure, and the num- 
ber of subsystems is not limited to two - and therefore no 
dual picture is present - the quest for a deterministic equation 
of motion for entanglement seems hopeless. For these cases a 



statistical approach is more encouraging [28]. 
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